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Abstract. In this paper we study synchrotron and synchrotron self Compton (SSC) emission from internal shocks (IS) during 
the prompt and X-ray flare phases of gamma-ray bursts (GRBs). The aim is to test the IS model for the flare emission and for 
whether GRBs can be GeV sources.We determine the parameters for which the IS model can account for the observed prompt 
and X-ray flares emission, and study the detectability of the high energy SSC emission by the AGILE and GLAST satellites.We 
find that the detectability of the SSC emission during the prompt phase of GRBs improves for higher values of the fireball 
Lorentz factor T and of the temporal variability t v . If IS is the mechanism responsible for the flare emission, and the Lorentz 
factor of the shells producing the flare is T ~ 100, the flare light curves are expected to present some substructures with temporal 
variability t v = 10 - 100 ms which are much smaller than the average duration of flares, and similar to those observed during 
the prompt phase of GRBs. If one assumes lower Lorentz factors, such as T ~ 10 -r 25, then a larger temporal variability t r ~ 40 
s can also account for the observed flare properties. However in this case we predict that X-ray flares do not have a counterpart 
at very high energies (MeV-GeV). An investigation on the substructures of the X-ray flare light curves, and simultaneous X-ray 
and high energy observations, will allow us to corroborate the hypothesis that late IS are responsible for the X-ray flares. 
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1. Introduction 

Observations with the Swift satellite have shown that the phase 
of prompt-to-afterglow transition is characterized by a very in- 
teresting phenomenology. Typically the tail of the prompt emis- 
sion, described by a very steep power law segment (temporal 
decay index 6 ~ 3 - 5), is followed hundreds to thousands of 
seconds later by a flat power law segment (6 ~ 0.5), and then 
around tens of thousands of seconds later by the usual after- 
glow power law emission (6 ~ 1.3). Sometimes there is a third 
break where the afterglow light curve steepen s to a power law 
6-2 dNousek et al. Il2006tlzhang et al. 1120061) . It has been ob- 
served that the tail of the prompt emission, as well as the flat 
power law phase and t he third segment are characterized by 



the p resence of flares dChincarini et al. I 120071: iFalcone et al. 



2007). Flares appearing hundreds of seconds after the burst 



were observed for the first time by BeppoSA X in a few events 



(e.g GRB011121, G RB011211, XRF01 1030. iPiro et al.ll2005 



Gall i & Piro 1 12006). Later Swift, thanks to its fast re-pointing 



capabilities, found that flares are a very common phenomenon 
in GRBs , as they are detect able in ~ 30 - 40% of its bursts 
sample ( O'Brien et al.l l2006). Moreover flare s are present both 



shift (e.g. GRB05 0803.lBloom et al. 112005b and high redshif t 
(e.g. GRB050904. [Cusumano et al. ll2007l:lGendre et al. 1120071) 
events. Several bursts show multiple flares in their X-ray light 
curves (e.g GRB 050607 and GRB 05073Q. |Pagani et al. 12 006: 



Perri et al. 1 120071) whil e others have o nly one big flare (see e.g. 
GRB060526 lchincarini et aTll2006l) . 



The flare phenomenon is complex, and despite the large 
number of ideas put forward in literature to explain its ori- 
gin, no model has completely interpreted the phenomenol- 
ogy of observed flares. The most important difference be- 
tween models is the duration of central engine activity: 
some of them require a long-lived cent ral engine, e.g. mod- 
els involving la te internal shocks (LIS) ( Burrows et al. 112 005: 
Wu et al. 112006b and delay e d exte rnal shocks (PES) dPiro et ail 



20051: iGalli&Pirol 120061 120071) produced by a long dura- 
tion and/or re-activation of the activity of the engine, while 
others do not require a long duration central engine activity, 
e.g. models involving LIS from a short duration central en- 

j 1 — i i 1 1 — — i* 

gine (Zhang et al. 2006; Wu et al. 2006) , refres hed shocks 
(Rees & Meszaros 1998; Kumar & Piran 2000a), two com- 



ponents jet (Meszaros & Ree s 1 120011 ; iLipunov et al. 1 1200 11) 



in lon g and short GRBs (e.g. GRB050724, iBarthelmv etaD P atch Y j et dKumar & Piran I l2000bl) forward shock (FS) 



2005b . in classical long G RB, X-ray rich (XRR ) and X-ray flash reverse shock (RS) dFan & Wei 1 120051). external shock (ES) 



-i^vi-a,, .„ „ p — ^ ^. y - J ^ II _| i , 1 

(XRF) (e.g. XRF050406 jRomano et al. 12006b . and in low red- with a clumpy medium dDermer & Mitotan | [1999j |Dermer 



Send offprint requests to: alessandra.galli@iasf-roma.inaf.it 



2007) and delayed magnetic dissipation in strongly magnetized 
ejecta caused by external shock dGiann ios 2006). 
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A temporal and spectral analysis of the first survey of 
X-ray flares observed by Swift has recently been carried out 



(Chi ncarini et al. I l2007t Falcone et al. 1 120071) . These authors 



analyzed the first 110 bursts observed by Swift and found that 
33 GRBs have significant flares in their X-ray light curves, 
for a total of ~ 70 flares. The temporal analysis revealed that 
there is a positive correlation between the time of flare ap- 
pearance and its duration, and that flares can be very ener- 
getic in compari s on wi th the underlying afterglow emission. 
Chincarini et al. 1 d2007l) find that a correlation exists between 
the ratio of intensity of successive prompt y-ray pulses and that 
between successive X-ray flares in the same burst. This is an 
indication that they m ay be caused by the same mechanism. 
Falcone et al. I (12007b used several spectral models to fit each 
flare. Usually flares can be fit ted both by a sim ple absorbed 
power law or by a Band model (Ba nd et al. Il993h . Some flares 
have a soft spectrum co nsistent with that of the af terglow emis- 
sion (e.g. GRB050712, lDe Pasquale et al. 1120061) . however the 
majority present a hard-to-soft spectral evolution. A spectral 
evolution of the order of 0.5-1.0 (depending on the index of 
the electron population) can be explained in the context of 
the ES if the typical emission frequency i s crossing the obser- 



vation al band, and models such as DES (Galli & Piro 2006, 



2007b . ES on a clumpy me dium (IDermer 



2007^ or refreshed 



shocks (IGuetta et al. 112007b can be applied to some flares of 
the sample. However a cco rding to the studies performed by 
Chincarini et al. I (120071) and lFalcone et al. I (12007b . a large frac- 



tion of flares could be explained by LIS produced by a long 
duration central engine. 

In this paper we focus on the possibility that flares are pro- 
duced by LIS and study the high energy emission from flares 
in this scenario. We plan to extend our study to other mod- 
els in a future study. In the LIS model the mechanism re- 
sponsible for the flare emission is the synchrotron radiation 
from relativistic electrons accelerated in the shocks. An ad- 
ditional radiation mechanism that may also play an important 
role is the synchrotron self Compton (SSC) emission, i.e. the 
up-scattering of the synchrotron photons by the same relativis- 
tic electrons to much higher energies. The synchrotron and SSC 
components from IS have been considered in the context of 



the prompt emission by IGuetta & Granot I (12003b . Following 
Guetta & Granot I (|2003|) . we estimate in this work the high en- 
ergy emission produced by SS C of the X-ra y flare p hotons. As 



noted by s everal authors e.g. Wang et al. I (12006b . iFan et al. 



d2007b and lGalli&Pirol (120071) . the detection of the predicted 
high energy (MeV to GeV) flare emission by observations with 
AGILE and GLAST combined with the X-ray flares detection 
by Swift would enable us to check the validity of different mod- 
els proposed to explain the afterglow phenomenology of GRBs. 

The paper is organized as follows. In Sect. 2 we report 
the AGILE and GLAST sensitivities. In Sect. 3 we summa- 
rize the calculations of the synchrotron and SSC spectra in 
Guetta & Granot (2003) for the prompt emission and provide 
the expression for the peak energy and the flux normalization. 
In Sect. 4 we show how the peak energies and flux normaliza- 
tion change for X-ray flares and estimate the high energy flare 
fluxes. Finally in Sect. 5 we give our summary and conclusions. 



2. Detector sensitivities 

In this section we present an estimate of the sensitivity of the 
high energy detectors aboard the AGILE and GLAST satel- 
lites. This is necessary to assess the detectability of high en- 
ergy emission during the prompt and flare phases of GRB we 
calculate in the following sections. 

The AGILE satellite is equipped with the GRID (Gamma- 
Ray Imaging Detector) instrument operating between 30 MeV 
and 50 GeV. We estimate the GRID sensitivity by adopting 
the criterion that a detection is made when at least 5 pho - 



tons are colle cted by the d e tector dZhang & Meszaros 1120011) . 
According to iGalli & Piro I (|2007) when the detector is source 
dominated, the detection threshold is given by: 



F, h (E) 



A eff (E)T b 



cm 2 s 1 



(1) 



where E is the photon energy, A e ff is the effective area of 
the detector at this energy, and r,„, is the detector integration 
time. For r„„=50 s the GRID is source dominated, thus we 
estimate the GRID sensitivity using the equation above. In cal- 
culating the GRID sensitivity we assume its effective area to be 
A e ff = 550 cm 2 throughout the energy range. 

The Large Area Telescope (LAT) aboard GLAST 
(Gamma-ray Large Area Telescope) observes between 
30 MeV and 300 GeV. The LAT effective area varies 
with energy in all the band; we take the value of the 
effective area at different energies from the LAT web page, 
http://www-glast.slac.stanford.edu/software/IS/glastJat_performance.htm 
As shown by IGalli & Piro (2007) for an integration time 
Ti„, - 50 s, the LAT is source dominated at all energies, thus 
we can estimate its sensitivity using Eq.Q] 



3. Synchrotron and SSC prompt emission from 
internal shocks 

In this section we estimate the synchrotron and SSC emission 
in the framew ork of the IS model, fo llowing the prescriptions 
presented in IGuetta & Granot I (12003b . In this model the flow 
Lorentz factor F is assumed to vary on a typical time scale t v 
and with an amplitude 6F ~ F. The shells collide at a radius 
R * 2F 2 cf, = 6 x 10 13 r| 5 ? v _ 2 cm, where F 2 . 5 = F/10 2 5 and 
?v,-2 = t v /(10~ 2 s). The internal energy released in each colli- 
sion is distributed among electrons, magnetic field and protons 
with fractions e e , e B and (1 - e e ) respectively. The electrons 
are accelerated in the shocks to a power law distribution of en- 



ergy N(y) cc y~P, and radiatively cool by the combination of 
synchrotron and SSC processes, the timescales of which are 
t syn ~ bnm e c I ctj B 2 y and t$sc = t syn /Y, the combined cooling 
time being t c = (l/t syn + l/tssc) 1 — + Y)> whe re B is the 



magn etic field, and Y is the Compton y-parameter dSari et al 
19961) . Y ~ e e /e B for e e « e B and Y « (e e /e B ) l/2 for e e » e B 
.The synchrotron spectrum is 
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where v sa , v ac and vm are defined in Guett a & Granot 



(2003) and the synchrotron peak energy is given by: 



E p = hv m = 0.12(1 +z)- 1 (l + Fr'^^Ljfr^f,;^ MeV 

(3) 

for p = 2.5. From the normalization condition we find the 
flux at the peak of the spectrum to be: 



Lis. 



24nD 2 



- = 1.3 X 10 _t, L 52 £> 2 8 erg cm s 



-2 „-l 



(4) 



where D 2 8 is the burst distance in unity of 10 28 cm. 
The SSC spectrum is given by 
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where Vj a c , v s ^ and v^, are defined in iGuetta & Granot I 



KN 



(2003) and the SSC peak energy is: 



E s p c = hv^ =4.6xl0 4 (l+z) 



^-L^K^MeV (6) 



for p = 2.5. 

For details about the spectrum at fr equencies higher 



than the Klein-Nishina frequency v^, see IGuetta & Granot 
(2003b). In our estimate of high energy emission we also take 
into acco unt the suppression of the emission due to pair pro- 
duction dGuetta & Granot Il2003h . 

In Fig. Q] we present the predicted synchrotron (dot-dashed 
lines) and SSC (dashed lines) prompt emission spectra as a 
function of the fireball Lorentz factor F for a burst with lu- 
minos ity L = 10 52 ere s located at redshift z=l. As in 
Guetta & Granot I d2003l) we take p=2.5, e e = 0.45, e B = 0.1, 



t v = 1 ms and vary the fireball Lorentz factor T between 200 
and 600. In the same figure we report the AGILE (solid line) 
and GLAST (dot-dot-dot-dashed line) sensitivities (see Sect. [2] 
for details on the AGILE and GLAST sensitivity). We note that 
th e spectral cutoff energy due to pair production (see Eq. (19) 
of IGuetta & Granot 1120031) moves to higher energies for higher 
fireball Lorentz factors while the peak of the prompt emission 
shifts to lower values. Therefore, given a set of parameters and 
a fixed integration time, the detectability of the SSC compo- 
nent improves with increasing fireball Lorentz factor values 
(see Fig. [T|i. The best candidate bursts for high energy detec- 
tion will be those with the peak energy close to the lower value 



of the BAT energy band, E p ~ 10 keV. We find that the high 
energy emission for a burst with these characteristics will be 
detected both by AGILE and GLAST if T > 350 for an integra- 
tion time Ti nt = 50 s and T ml = 10 s. In the most favorable case, 
F ~ 600 and T,„, = 50 s , the SSC prompt emission component 
can be detected up to maximum redshift z mfll ~ 3.2 by AGILE, 
andz m ax ~5.5 by GLAST. 



L= 1 52 erg s~\ €,=0.45, e„=0.1, p = 2.5, z=1.0 
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Fig. 1. Synchrotron (dot-dashed lines) and SSC (dashed lines) 
prompt emission spectra for a burst at redshift z = 1. The 
burst prompt emission luminosity is £52=1 and its variability is 
ty = 1 ms. We fix the fraction of fireball energy going into rel- 
ativistic electrons and magnetic field to e e = 0.45 and eg = 0.1 
respectively, and fix p — 2.5. We display the predicted spectra 
for three different values of the fireball Lorentz factor F: 200 
(purple), 350 (blue) and 650 (red). The solid lines and the dot- 
dot-dot-dashed lines represent AGILE and GLAST sensitivity 
for an integration time of 10 s (in grey) and 50 s (in black). The 
solid vertical lines refer to the Swift XRT (red) and BAT (black) 
energy ranges. 



As in IGuetta & Granot I d2003l) we also study the predicted 
synchrotron and SSC emission as a function of the burst tem- 
poral variability t v for a burst at z = 1.0 (Fig. [2]). In this case, 
we fix F = 600. The prompt SSC peak energy (Eq.|6]l moves to 
lower energy with larger t v while the spectral energy cutoff due 
to pair production (Eq. 19 of lGuetta & Granot 2003) moves to 
higher energies with larger t v , therefore the detectability of the 
SSC prompt emission component improves with an increasing 
variability time f v . 

Figure|2]shows that the SSC component can be detected by 
AGILE and GLAST for t v = 0.1-10 ms both for 10 s and 
50 s of integration time, and that the best condition for the de- 
tectability of the SSC component is achieved for t v = 10 ms. 
Therefore also in this case the best candidate bursts for high 
energy detection would be those with E p close to the lower 
value of the BAT energy band. For t v = 10 ms and an instru- 
mental integration time r,„, = 50 s, the SSC prompt emission 
component can be detected up to maximum redshift z max ~ 6.5 
by AGILE, and up to z ma x ~ 7.0 by GLAST. We summarize 
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our results in Table [TJ where we report the maximum redshift 
at which AGILE and GLAST can detect a burst as a function 
of the fireball Lorentz factor T, the prompt emission temporal 
variability t v and the instrumental integration time T im . Larger 
values of T and t v shift the cutoff due to pair production to 
higher energies and the peak energy of the prompt emission 
to lower energy. Actually the peak energy of the synchrotron 
emission of these GeV emitters bursts falls below or inside the 
(15-150 keV) BAT band of Swift (see Fig.Q). For this purpose, 
the energy band covered by Swift is better suited than the band 
covered by BATSE (300-500 keV) for the identification of the 
GeV burst sources. 



Table 1. Maximum redshift of detection by AGILE and 
GLAST as a function of the fireball Lorentz factor F, the 
prompt temporal variability t v , and the instrumental integra- 
tion time Tint. The other parameters are fixed as L52 = 1.0, 
e e = 0.45, €b = 0.1, and p = 2.5. We also report the peak 
energies E p of the synchrotron prompt emission and the cutoff 
energies E cut of the relative SSC component. 



10" 5 



10"" r 



10" 7 . 



L= 1 52 erg s~\ e e =0.45, e B =0.1 , p = 2.5, z=1.0 



t„=0.1 ms t„=1 ms t„=10 
r=600 




AGILE 



10" 4 



10" 



10° 1 o 2 
hi/ [MeV] 



10" 



10 s 



Fig. 2. Synchrotron (dot-dashed lines) and SSC (dashed lines) 
prompt emission spectra for a burst at redshift z = 1.0. The 
fireball Lorentz factor is fixed to T — 600, and the other param- 
eters are the same as in Fig.Q] We display the predicted spectra 
for three different values of the burst temporal variability t v : 
0.1 ms (red), 1.0 ms (blue) and 10 ms (purple). The solid lines 
and the dot-dot-dot-dashed lines represent AGILE and GLAST 
sensitivity for an integration time of 10 s (in grey) and 50 s (in 
black). The solid vertical lines refer to the Swift XRT (red) and 
BAT (black) energy ranges. 



4. X-ray and high energy flares in the context of 
internal shock 

As stated in Sect. Q] several models have been proposed in 
the literature to explain the origin of flares. Some GRBs 
show properties which can originate from LIS produced b 
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sion radii and consequently of all emission properties of X- 
ray flares. All the relevant quantities, such as the synchrotron 
(Eq. [3]l and SSC (Eq. |6]) peak energy, with the correspond- 
ing fluxes and the high energy cutoff, depend strongly on t v . 
If the mechanism responsible of the X-ray flares is the same 
as the prompt emission (IS) one could expect a similar vari- 
ability time, thus we initially assume for flares that t v ~ 100 
ms. However, given our ignorance on the average interval be- 
tween consecutive shells when the engine is reactivated, that 
these shells could be emitted close to the central engine or at 
larger distances, and that the Lorentz factor F could be smaller 
for X-ray flares in comparison with the prompt emission, in the 
following we also discuss the case of a larger t v . The other rele- 
vant quantities are the luminosity L of the flare and the average 
flare duration tf. In order to determine the typical tf and L of 
the flares, we select a sample of X-ray flares with known red- 
shift and published light curves, taking place at times shorter 
than » 1000 s, when the IS mechanism can still be active. 
At longer times other mechanisms can be responsible for the 
flare emission, e.g. inh omogeneities in the external medium 
dLazzati & Perna ll2007b . 



We take our sample of flares from Chi ncarini et al. 1 (12007) 



a lon g lasting centr al engine activity (IChincarini et al. 11200 



Falcone et al 



2007). This mot i vated us to apply the IS model 



presented by lGuetta & Granot I (|2003l) to the flare phenomenon 
and study the GeV emission from flares. In order to estimate 
the high energy counterpart of X-ray flares, we repeat the cal- 
culations given in the previous section. 

As shown in Sect. [TJ in the LIS scenario the variability 
time t v (i.e. the average interval between consecutive shell ejec- 
tions) has a fundamental role in the determination of the colli- 



in Table [2] w e report the flare prope rties which are relevant 
to our study. IChincarini et al. (2007) performed the temporal 
analysis of a sample of 33 GRBs detected by Swift using a 
multi-broken power law to fit the underlying continuum and 
a Gaussian model for X-ray flares, and adopted the Gaussian 
width cr as a measure of the X-ray flare duration tf. In order to 
estimate the high energy emission which can be associated to 
the X-ray flares, we need the flare peak luminosity which we 
can get from the temporal and spectral information available in 
literature. 
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In Fig. [3] |U and [6] we report the flare peak time, du- 
ration, redshift and luminosity distributions of our sample of 
flares. For comparison we also report in the same figures the 
distributions of the complete (with a nd without known red - 
shift) sample of flares presented by IChincarini et al. I (2007). 



As we can see from these figures, the total sample of flare 
and our sub-sample of bursts with known redshift show the 
same distributions. For the total sample, the mean peak time 
of flares is (t p ) = 404.0 s, cr = 300.5 and the mean duration is 
(tf) = 40.4s, cr = 41.8 s. We get similar values for the bursts 
of our sam ple, (t„) = 461.5 s, cr = 330.6, and (tf) = 39.6s, 
cr = 27.6 dChincarini et al. 112007b (note that the fit of flares 



with two power laws gi ves flare durations o f the order of hun- 
dreds of seconds, e.g. Falcone et al. 1 12007). We take the peak 



flux of each flare from published light curves and find their 
logarithmic mean peak luminosity, (logL p ) = 49.3, cr = 1.2 i.e 
L p ~ 2x 10 49 erg s~ l , while if we compute the "linear" medium 
we find (U) ~ 10 50 erg s _1 . The mean redshift of the bursts of 
our sample is (z) = 2.76, cr - 1.92. 



Peak time X- 



20 



15 



10 



■ray 



flares distribution 




600 800 1000 1200 1400 
t p [s] 

Fig. 3. Distribution of peak times of the total (with and with- 
out redshift) flare sample of IChincarini et alT (2007) (in grey), 
and of the sub-sample of flares with known redshift (in purple) 
reported in Table [2] 



To estimate the high energy counterpart of X-ray flares, we 
repeat the calculations made in Sect. Q]by changing only the 
luminosity fro m L = 10 52 ergs~ l to L = 2 x 10 49 



erg s 



As noted by lFan et al. I (12007) it is difficult to predict the ex- 
pected SSC emission as we do not have a good estimate of the 
typical Lorentz factor of the electrons accelerated at the shock 
to produce the X-rays. However we can give some constraints 
on the Lorentz factor T of the colliding shells. If we assume that 
the X-ray flares are produced by accelerated electrons which 
cool very fast, the condition v c < v m introduces an upper limit 
to the fireball Lorentz factor F. In particular, for typical flare lu- 
minosity and temporal variability f,,=100 ms, and for e e = 0.45, 
£b = 0.1 and p = 2.5, we find F < 340. In addition, in order to 
have SSC GeV flare emission we have to require that the spec- 
tral cutoff energy due to pair production is > 1 GeV. With the 
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30 
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300 



Fig. 4. Distrib ution of the flare durati on of the complete sample 
of flare from lChincarini et al. I d2007l) (in grey), and of the sub- 
sample of flares with known redshift reported in Table [2] (in 
purple). 



redshift distribution 




Redshift 



Fig. 5. Redshift distribution of the flares of our sample. 



above quantities, using Eq. 19 of iGuetta & Granot I Q2003), we 
find T > 60. We thus take in the following r = 100 and T = 300 
as two possible Lorentz factor values. The peak energy of syn- 
chrotron and SSC flare emission are given by: 



Epjiare = 3.8(1 + z^l + Y)~ 1/3 ^e^L^T^t^ keV (7) 



and 



7 sc 

J p, flare 



= 1.5(1 + zr l e:' 2 e l J 2 L l 4 i 2 r- 2 2 C l GeV 



(8) 



with L49 the flare peak luminosity in unity of 10 49 erg s _1 , 
T2 the flow Lorentz factor in unity of 100, and t v -\ the flare 
temporal variability in unity of 100 ms. For typical parame- 
ters values the peak of flare synchrotron emission is inside or 
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Table 2. Our sample of X-ray flares with known redshift and published light curves. t p is the time of the X-ray flare peak, tf 
is the flare duration, At/t is the duration-to-flare time ratio, and L p is the peak l uminosity of X-ray fla res e valuated as L n - 
J ' - kw(1 + z) (2 ~ y) wit h y t he power law pho ton i ndex. References: [ 11 iRomano et al. I (12006): [2] IChincarini et al. 



4nDfF .2- 



1 ).o , 



(2007); [3] Falcone et al. I (120071) ; 141 iPerri et al. I (120071) ; r51 ICenkoet al. 
(2006). 



(2006) ; [61 ICusumano et al. I d2007): [71 lQates etal. 



name 


redshift 




f f 


At/t 




References 




z 


[s] 


[s] 




[erg s ] 




GRB050406 


2.44 


211 ±5 


17 9 +12J 


0.882 


2.3 x 10 49 


[1] 


GRB050724 


0.26 


275 ± 5 


30 6 +6 6 




3.5 x 10 47 


T21 HI 






327 + jJ 


12 7 + " 

' -5.0 




1.6 x 10 47 




GRB050730 


3.97 


131.8™* 
234.2+2-7 

436.5^ 


14.5+3 g 
32.7+fj 4 




1.0 x 10 3U 
8.3 x 10 49 

1.1 x 10 50 


[4] 






685.8!^ 






7.8 x 10 49 




GRB050802 


1.71 


464 + 31 


100 +33 


2.327 


4.3 x 10 48 


[2],[3] 


GRB050820A 


2.61 


241« 


9 5 +U3 




1.2 x 10 5U 


[5] 


GRB050904 


6.3 


448.6«- 7 


45.9 + ^ 




8.1 x 10 5U 


[6] 




6.3 
6.3 


975 - 5+ 32'5 

1265.5+^ 


62.8^6-9 




8.5 x 10 49 
8.7 x 10 49 




GRB060108 


2.03 


303.5^j 


44 c+ 125.5 
-30.5 


1.405 


4.0 x 10 47 


[7] 



X-ray flares peak luminosity distribution 




L — 2 x 10 49 erg s~' the (predicted) high energy flare emission 
can be detected by AGILE and GLAST up to a maximum red- 
shift Zmax ~ 1 .2 and z max ~ 1 .4 for an integration time of 500 s, 
and up to Zmm, ~ 0.6 and z ma x ~ 0.7 for an integration time of 
100 s. 



L=2*10 49 erg s~\ £.=0.45, e B =0. 1 , p = 2.5, z=1.0 



0.0001 0.0010 0.0100 0.1000 1.0000 10.0000 100.0000 
Peak Luminosity L p [10 erg s~'] 

Fig. 6. Luminosity distribution of the flare of our sample. 



just below the XRT band and the peak of flare SSC emission is 
around ten-hundred MeV. The shells producing the flares col- 
lide at a radius R ~ 6 x 10 13 cm from the central engine for 
T = 100, i.e close to the typical distances where the collisions 
producing the prompt emissio n take place. This is an important 



r=ioo r=300 

t =100 ms 



SSC 




/ LAT 



10" 



10"* 



10" 2 10° 10 2 

hi/ [MeV] 



10" 



10" 



result because, as stressed by iFan et al. I (120071) the knowledge 



of the location of the shocks determines the parameters of the 
emission, which strongly affect the estimation of the SSC com- 
ponent. 

As for the prompt emission we study the predicted syn- 
chrotron and SSC flare emission as a function of the fireball 
Lorentz factor T, at redshift z — 1 (Fig. |7j>. The flare duration is 
estimated to vary from tens to hundreds of seconds depending 
on the model used for the fit (i.e. a Gaussian or a power law 
respectively), we thus assume for AGILE and GLAST integra- 
tion times of 100 s and 500 s. We find that for a flare luminosity 



Fig. 7. Synchrotron (dot-dashed lines) and SSC (dashed lines) 
spectra for an X-ray flare with luminosity L = 2x 10 49 erg s~ l 
and temporal variability t v - 100 ms at redshift z - 1 as a 
function of the fireball Lorentz factor T. Other model parame- 
ters are the same as in Fig.Q] We display the predicted spectra 
for T=100 (red) and T=300 (blue). The solid and dot-dot-dot- 
dashed lines represent AGILE and GLAST sensitivity for an 
integration time of 100 s (in grey) and 500 s (in black). The 
solid vertical lines refer to the Swift XRT (red) and BAT (black) 
energy ranges. 
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We then calculate the emission expected for X-ray flare lu- 
minosity L' = 10 50 erg s _1 as a function of the fireball Lorentz 
factor r, at z = 1 (Fig. [8j. With this X-ray flares luminosity 
the predicted high energy emission can be detected both by 
AGILE and GLAST at z = 1 . If we take as mean flare luminos- 
ity (L') = 10 50 erg s -1 , the SSC flare emission can be detected 
by AGILE and GLAST up to maximum redshift z ma x ~ 2.4 and 
Zmax ~ 2.8 for an integration time of 500 s, and up to Zmax ~ L2 
and Zmax ~ 1.4 for an integration time of 100 s. 
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Fig. 8. Synchrotron (dot-dashed lines) and SSC (dashed lines) 
spectra for an X-ray flare with luminosity L = 10 50 erg s~ l and 
temporal variability t v = 100 ms at redshift z — 1 as a func- 
tion of the fireball Lorentz factor T. The other model param- 
eters are the same of Fig. Q] We display the predicted spectra 
for T=100 (red) and F=300 (blue). The solid and dot-dot-dot- 
dashed lines represent AGILE and GLAST sensitivity for an 
integration time of 100 s (in grey) and 500 s (in black). The 
solid vertical lines refer to the Swift XRT (red) and BAT (black) 
energy ranges. 



We also study how the relative importance of synchrotron 
and SSC emission varies with the ratio (e e /es), this ratio be- 
ing the quantity which determines the relative importance of 
these two emission mechanisms (for a detailed analysis of the 
effects of this ratio during early and late afterglow emission see 
iGalli & Piro 1120071) . We present our results in Fig. |9]for X-ray 
flare luminosity L — 2 x 10 49 erg s' 1 . As expected, the relative 
importance of the two processes increases with the ratio (e e /eg) 
thus favoring the detection of the high energy component even 
though the peak of SSC emission moves toward lower energies 
with smaller eg values. We keep e e = 0.54 and vary eg from 0.1 
(red curves), to 0.01 (blue curves), and finally to 0.001 (pur- 
ple curves). The predicted SSC emission can be detected by 
AGILE and GLAST up to z max ~ 2.5 and z ma x ~ 2.7 for an inte- 
gration time of 500 s, and up to z max ~ 1 .2 and Zmax ~ 1.3 for an 
integration time of 500 s. For a luminosity L = 10 50 erg s~ l the 
high energy emission can be detected by AGILE and GLAST 
up to z,nax ~ 5.2 and z ma x ~ 5.5 for an integration time of 500 



s, and up to Zmax ~ 2.6 and Zmax ~ 2.8 for an integration time 
of 100 s. 
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Fig. 9. Synchrotron (dot-dashed lines) and SSC (dashed) spec- 
tra for an X-ray flare with luminosity L = 2 x 10 49 erg s~' and 
temporal variability t v - 100 ms at redshift z — 1 as a function 
of eg- The fireball Lorentz factor is fixed to T — 100, and other 
model parameters are the same as in Fig.Q] We display the pre- 
dicted spectra for 6b=0.1 (red), e#=0.01 (blue) and e#=0.001 
(purple). The solid and dot-dot-dot-dashed lines represent the 
AGILE and GLAST sensitivities for an integration time of 100 
s (in grey) and 500 s (in black). The solid vertical lines refer to 
the Swift XRT (red) and BAT (black) energy ranges. 



Falcone et al. I (120071) have shown that X-ray flares spectra 
can be fitted both by a power law and/or by a Band model. 
However, due to the incomplete spectral coverage they could 
constrain the flare peak energy only in some cases, and sug- 
gested that typically this peak energy should be in the soft X- 
ray band, or between the optical and X-ray bands. Figures [7] 
[8] and [9] show that if we assume t v ~ 100 ms and T ~ 100, 
the peak energy of the X-ray flare emission falls in, or just 
below, the XRT band consisten t with the findings of the flare 



spectral analysis performed by iFalcone et al. I (120071) . In Fig. 



JU]we show the predicted synchrotron and SSC flare emission 
for a range of t v values similar to that typically observed dur- 
ing the prompt emission, i.e. t v ~ 10 ms - 1 s, with T = 100 
and L — 2 x 10 49 erg cm" 2 . As we can see from this figure, 
when one assumes a Lorentz factor of the order of one hun- 
dred for this range of t v the peak energy of the X-ray flare is 
always inside or just below the XRT band. The detection of the 
SSC component related to the X-ray flare depends strongly on 
the flare temporal variability: for t v < 10 ms the spectral cutoff 
due to pair production is shifted to lower energies (Eq. 19 of 



Guet ta &Granot 1 12003). i.e. in the AGILE and GLAST bands, 



making the SSC emission more difficult to detect. In the same 
way, a larger temporal variability t v > 1 s shifts the peak energy 
of the SSC component below the AGILE and GLAST energy 
bands, and this again makes the high energy flare more difficult 
to detect. In order to have both the X-ray peak energy in the 
XRT band and a detectable high energy emission for F ~100, 
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we need ?,, ~ 50-100 ms. This implies that the flux of the X-ray 
flares should vary on a time scale much smaller than the flare 
duration, i.e t v « tj. Therefore if LIS is the mechanism re- 
sponsible for the flare emission, and the flare has a counterpart 
at high energies (from MeV to GeV), the X-ray flare light curve 
should present some substructures where the flux can even dou- 
ble its value at a time t v which is much smaller than the flare 
duration as observed in the prompt emission. 

We should note however that the temporal resolution of 
the Swift XRT is ~ 2 ms in WT mode and 2.6 s in PC mode 
(see http://www.swift.psu.edU/xrt/software.html#modes). At 
the time of flare occurrence the readout is typically in the PC 
mode. Data are taken in WT mode only for the brightest flares, 
however the search of substructures with duration < 1 s would 
be important to check the LIS model. If no significant varia- 
tion in the X-ray flux on a small time scale is detected, then 
the average time interval t v between consecutive shell ejections 
should be the average value of the flare duration. For our sam- 
ple we obtain f v ~ 40 s, which would change the lower and 
upper limits of the flow Lorentz factor F to ~ 25 and ~ 110 
respectively. 

If one assumes the Lorentz factor to be of the order of 50, 
the LIS model predicts that the peak energy of the X-ray flare 
emission is around 1 eV, i.e. well below the XRT band (see Fig. 
Hfl. For t v = t f = 40 s, L = 2 X 10 49 erg cm- 2 and F=50, the 
shells producing the flare would collide at a radius R ~ 6 x 10 15 
cm from the central engine, i.e. at distances larger than that 
where the collisions producing the prompt emission take place. 
For f,,=40 s if one assumes a lower Lorentz factor as F=25 or 
T =10, then the X-ray flare peak energy moves to higher values 
consistent with those found for t v ~ 100 ms. Under these as- 
sumptions large temporal variability values could also account 
for the observed flare properties. Moreover, low values of the 
Lorentz factor as T ~ 10 have been used by several authors, e.g. 
Falcon e et al. (2006), to fit X-ray flares. However in the case 



L=2.0*10 49 erg s~\ e,=0.45, e,=0.1, p = 2.5, z=1.0 



of low Lorentz factor values the spectral cutoff due to pair pro- 
duction goes below 1 GeV, and X-ray flares are not expected 
to have a high energy counterpart. Simultaneous observations 
by Swift in the X-ray band, and by AGILE and GLAST at high 
energies, and a detailed temporal analysis of flares light curves 
are thus very important to test the LIS model we propose in this 
paper. 

We summarize our redshift study for SSC flare emission in 
Table [3] where we give the maximum redshift achievable by 
AGILE and GLAST for an integration time T- mt — 100 s and 
Ti„t = 500 s as a function of the fireball Lorentz factor F, the 
flare luminosity L, the ratio (e c /eg), and the flare temporal vari- 
ability t v . The peak energy of flares moves to higher energies 
with increasing flare peak luminosities and lower values of the 
(e c /efi) ratio, while it shifts to lower energies for larger flow 
Lorentz factors. 

As can be seen from Table [3] if LIS are responsible for X- 
ray flares and flare light curves contain some features where the 
flux can double its value at f„ = 10-100 ms, similar to which is 
observed during the prompt emission, the peak energy of X-ray 
flares is in, or just below, the XRT e nergy band consistent with 
the spectral analysis performed by Falcone et al. (2007). For 
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Fig. 10. Synchrotron (dot-dashed lines) and SSC (dashed) spec- 
tra for an X-ray flare with luminosity L — 2 x 10 49 erg s as a 
function of the temporal variability t v at redshift z = 1 . The fire- 
ball Lorentz factor is fixed to T = 100 and other model param- 
eters are the same as in Fig.Q] We display the predicted spec- 
tra for f v =10 ms (red), f,,=100 ms (blue) and t v =l s (purple). 
The solid and dot-dot-dot-dashed lines represent the AGILE 
and GLAST sensitivities for an integration time of 100 s (in 
grey) and 500 s (in black). The solid vertical lines refer to the 
Swift XRT (red) and BAT (black) energy ranges. 
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Fig. 11. Synchrotron (dot-dashed lines) and SSC (dashed lines) 
spectra for a typical X-ray flare with luminosity L — 2 x 
10 49 erg s~ l and temporal variability t v = 40 s at redshift 
Z = 1 as a function of the fireball Lorentz factor T. Other 
model parameters are the same as in Fig.Q] We display the pre- 
dicted spectra for F=10 (red), T=25 (blue) and F=50 (purple). 
The solid and dot-dot-dot-dashed lines represent AGILE and 
GLAST sensitivity for an integration time of 100 s (in grey) 
and 500 s (in black). The solid vertical lines refer to the Swift 
XRT (red) and BAT (black) energy ranges. 



t v — 40 s the peak energy is expected to be between the optical 
and the X-ray band only for low Lorentz factor values. 
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Table 3. Maximum redshift of detection by AGILE and GLAST as a function of the fireball Lorentz factor F, the flare luminosity 
L, the ratio (fe/eg), and the flare temporal variability f v , and the instrumental integration time Ti nt . 



satellite 


r 


L 


e e le B 


ty 




<-max 


Ep 




Ecu; 






[ergs- 1 ] 




[ms] 


[s] 




[keV] 




[MeV] 


AGILE 


100 


2 x 10 49 


4.5 


100 


500 


1.2 


0.16 




8.0 x 10 3 


GLAST 


100 


2 x 10 49 


4.5 


100 


500 


1.4 


0.15 




7.4 x 10 3 


AGILE 


100 


2 x 10 49 


4.5 


100 


100 


0.6 


0.22 




1.1 x 10 4 


GLAST 


100 


2 x 10 49 


4.5 


100 


100 


0.7 


0.22 




10 4 


AGILE 


100 


10 50 


4.5 


100 


500 


2.4 


0.23 




1.2 x 10 3 


GLAST 


100 


10 iu 


4.5 


100 


500 


2.8 


0.21 




1.1 x 10 3 


AGILE 


100 


10 iu 


4.5 


100 


100 


1.2 


0.36 




1.9 x 10 3 


GLAST 


100 


10™ 


4.5 


100 


100 


1.4 


0.33 




1.7 x 10 3 


AGILE 


100 


2 x 10 49 


450 


100 


500 


2.5 


5.2 x 10- 


3 


8.0 x 10 3 


GLAST 


100 


2 x 10 49 


450 


100 


500 


2.7 


4.9 x 10- 


3 


7.6 x 10 3 


AGILE 


100 


2 x 10 49 


450 


100 


100 


1.2 


8.3 x 10- 


5 


1.3 x 10 4 


GLAST 


100 


2 x 10 49 


450 


100 


100 


1.3 


8.0 x 10- 


3 


1.2 x 10 4 


AGILE 


100 


10 iu 


450 


100 


500 


5.2 


6.6 x 10- 


3 


1.1 x 10 3 


GLAST 


100 


10 50 


450 


100 


500 


5.5 


6.3 x 10- 


3 


10 3 


AGILE 


100 


10 50 


450 


100 


100 


2.6 


1.1 x 10- 


2 


1.8 x 10 3 


GLAST 


100 


10 3U 


450 


100 


100 


2.8 


1.1 x 10- 


2 


1.7 x 10 3 


AGILE 


25 


2 x 10 49 


4.5 


4x 10 4 


500 


0.8 


7.8 x 10 


-3 


1.7 x 10 3 


GLAST 


25 


2 x 10 49 


4.5 


4x 10 4 


500 


0.9 


7.4 x 10 


-3 


1.6 x 10 3 


AGILE 


25 


2 x 10 4y 


4.5 


4x 10" 


100 


0.35 


0.01 




2.2 x 10 3 


GLAST 


25 


2 x 10 49 


4.5 


4x 10 4 


100 


0.4 


0.01 




2.2 x 10 3 



5. Summary and conclusions 

We have calculated the synchrotron and SSC emission during 
the prompt and the X-ray flare GRB phases from IS and studied 
how they depend on model parameters. For p — 2.5 (the power 
law index of the electron energy distribution) and typical model 
parameters, the SSC component dominates the prompt emis- 
sion above » 100 MeV and th e flare emission above » 100 keV. 
As can be seen from Eq. 1 9 o flGuetta & GranoT ( 2003 ) and Eq. 
[3] and [6] and from Fig. Q] and El during the prompt phase larger 
values of the flow Lorentz factor F and the temporal variability 
t v shift the cutoff energy to higher values, and the peak of the 
emission to lower values. For example, in order to have ~ 1 
GeV emission for t v - 1 ms, we need F > 350 which implies 
E p < 100 keV exactly in the range of BAT (Swift). Therefore 
we expect that a good fraction of the bursts detected by Swift 
can produce GeV emission which will be detected by AGILE 
and GLAST. 

We have taken the LIS as the mechanism responsible for the 
flare emission. In this model the average interval between con- 
secutive shell ejection f v , has a fundamental role in determining 
the radius of the collision and therefore the emission properties 
(i.e. the peak energy) of X-ray flares. We initially considered 
a range of t v (10 ms-1 s) and a Lorentz factor (r ~100) sim- 
ilar to the prompt emission, and found that the peak of the 
emission is between the optical and the X-ray bands, consis- 
tent with the observations. This implies that the X-ray flare 
light curve should contain some substructures of t v ~ 100 
ms such as the prompt emission. An accurate search of these 
flux variations is mandatory to understand the physics of the 
flares and to check the LIS model. If no flux variation on a 
small time scale is detected the implication would be that the 



shells are ejected with an average interval similar to the dura- 
tion of the single flare (?/ ~ 40s). In this case the collisions 
betwee n the shells happen at a large distance from the central 
engine. lSpada et al. I (120001) have shown that the pulse duration 



increases with the distance R from the central engine, and con- 
sequently with time. This is consistent with the observations 
which show that t he flare duration incr eases with the time of 
flare appearance (IChincarini et al. 112007b . If t v — tf ~ 40 s the 
peak energy of the X-ray flare emission is between the optical 
and the X-ray band only for low Lorentz factors, F ~ 10 -r 25. 
However, in this case the cutoff due to pair production is below 
1 GeV and X-ray flares are expected to not have a counterpart 
at high energies. 

We have studied the detectability of the synchrotron and 
SSC emission during the prompt and the X-ray flare GRB 
phases as a function of the model parameters. During the GRB 
prompt phase we find that the detectability of the SSC emis- 
sion component improves with increasing fireball Lorentz fac- 
tor values T and increasing t v . For a typical prompt luminosity 
L = 10 52 erg s _1 the SSC prompt emission component can 
be detected by AGILE and GLAST up to a maximum redshift 
Zmax ~ 3.0 and z„ W x ~ 4.0 for an integration time T int — 10 s, 
and up to z nmx ~ 6.5 and z m ax ~ 7.0 for an integration time 
T int = 50 s. 

During the X-ray flare phase the detectability of the SSC 
component improves with smaller values of the fireball Lorentz 
factor F depending on the position of the SSC peak energy and 
the cutoff energy due to pair production with respect to the ob- 
servational band of the detector, and with larger values of the 
fe/ee) ratio. For a flare luminosity L — 2 x 10 49 erg s and 
t v = 100 ms the SSC component can be detected by AGILE and 
GLAST up to a maximum redshift z m ax ~ 2.5 and z max ~ 2.7 
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for an integration time of 500 s, and up to z m ax ~ 1.2 and 
Zmax ~ 1.3 for an integration time of 100 s. If we assume a 
more optimistic mean flare luminosity L — 2 x 10 50 erg s _1 the 
SSC component can be detected by AGILE and GLAST up to a 
maximum redshift z max ~ 5.2 and z max ~ 5.5 for an integration 
time of 500 s, and up to and up to z, mx ~ 2.6 and Zmax ~ 2.8 for 
an integration time of 100 s. 

Our findings show that AGILE and GLAST have good pos- 
sibilities for detecting the high energy emission coming from 
flares. In particular, a comparative study of X-ray synchrotron 
emission and high energy SSC emission is an important tool 
to test the X-ray flare model we present in this paper, and/or 
to disentangle between the several models proposed in the lit- 
erature to explain the origin of flares. In the context of LIS 
models, an additional mechanism producing high energy flares 
is possible (with respect to the SSC process considered in this 
paper), i.e the external inverse Compton scattering of the X-ray 
flare photons on the afterglow electrons. We can discriminate 
this component from the SSC, because in this case we expect 
that the high energy flares last much longer than the X-ray flare 
dWang et al. 1 12006; Fan et al. 112007b . A distinctive element be- 



tween the variant of the IS model which we study in detail in 
this paper and mod els in which flares are caused by an ES (e.g 
Gall i & Piroll2007l) is that in the context of LIS the peak of the 
SSC flare component is expected to be at lower energies with 
respect to that of the ES SSC component, i.e. in the MeV band 
for IS and in the GeV-TeV band for ES. 

Acknowledgements. The authors are grateful to L. Stella for useful 
comments and discussions. 
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